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ABSTRACT 

Taking advantage of the smoothness of CMB Ce power spectra, we derive a simple 
and model- independent parameterization of their measurement. It allows to describe 
completely the spectrum, i.e. provide an estimate of the value and the error for any 
real / point at the percent level, down to low £ multipole. We provide this parameter- 
ization for WMAP first year data and show that the spectrum is consistent with the 
smoothness hypothesis. We also show how such a parameterization allows to retrieve 
the Ce spectra from the measurement of Fourier rings on the sky (Tm) or from the 
angular correlation function {C{9)). 
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1 INTRODUCTION 

The study of the anisotropies of the cosmic microwave back- 
ground (CMB) temperature has become a major field in 
modern cosmology. For angular scales below < 1°, the 
Ce power spectrum, defined as the Legendre-transform of 
the two-point auto-correlation function C{6), is expected to 
present a set of acoustic peaks due to causal physics well es- 
tablished on the theoretical ground {e.g. Hu and Dodelson 
2002). 

When the first peak began to emerge from the data, it 
was natural to characterize its location/amplitude/spread. 
This was performed by estimating the maximum in a fixed 
I range, for instance by fitting a gaussian (Knox and Page 
2000) or a polynomial function (Durrer et al. 2003). It was 
further motivated by the fact that Hu et al. (2001) showed 
that for the first two peaks most of the cosmological model 
information was contained in the peaks locations and rela- 
tive heights. 

As more peaks became available, in particular in the 
high I region thanks to interferometer-based experiments, 
a more complete phenomenological fit was proposed to de- 
termine the peaks locations (Odman et al. 2003) through a 
sum of gaussian functions. Adding an oscillatory function 
has also been proposed (Douspis and Ferreira 2002) to char- 
acterize the existence of the peaks. While acceptable in the 
50 < ^ < 1250 region (with 5 gaussians) it fails to describe 
the low i Sachs- Wolfe plateau and the high part of the spec- 
trum (Odman 2003). Note that there are no physical reasons 
for the peaks to have a gaussian shape. 

With the advent of the high precision WMAP results 
(Bennet et al. 2003) and the expected huge sensitivity of 



the future planned Planck satellite ^ mission, it is time to 
consider the precise parameterization of the full Ce spectrum 
over a broad £ range. 

Obviously a modeling through cosmological parameters 
is not adapted to describe a purely experimental spectrum. 
We propose to take advantage of the expected smoothness 
of the Ce power spectrum, which comes from a combined 
effect of the continuity of the Fourier P{k) spectrum in the 
Standard Model, and the use of spherical Bessel functions 
to project it onto the £ space {e.g. Bartlett 1999). 

This smoothness property has already been exploited 
for fitting the Ce spectrum with splines {e.g. Oh et al. 1999). 
Here we will rather work in the light of the Fourier decom- 
position in particular by revisiting the sampling theorem. 
This will allow to provide a very simple description of any 
Ce spectrum as a function of real £ values. Therefore our 
goal is twofold: 

(i) obtain a reduced number of parameters to describe a 
CMB power spectrum 

(ii) provide an interpolation for any (real) £ value, i.e. a 
central value and a (Gaussian) error. 



2 DESCRIBING Ce SPECTRA 

The sampling theorem A CMB power spectrum Ce = 
£{£+ l)/2n Ce can be considered as the sampling for integer 
values of a smooth continuous signal C{1), where I is real. 
The signal is band-limited so that its Fourier transform can 
be neglected for |/| > /max. The sampling theorem states 
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Figure 1. Three different CMB TT power spectra (up to ^max = 
3000), used as illustrations for the method. The darlc full line 
corresponds to WMAP data best fit model. 



that for any rate (A^) above the critical frequency: 
1 



the real space complete signal can be recovered through: 



(1) 



e{l) = ^ Cj sine ( TT 



(2) 



where Ij are the sample positions. 

It appears that "reasonable" CMB power spectra are 
(mostly) band-limited to low frequencies. We illustrate that 
feature by choosing three representative spectra (Figure 1): 
one is the "WMAP best fit model" {AC DM of Spergel et 
al. (2003)), a second one has different peak proportions, and 
a third one has shifted peak locations. 

We then show their power spectral densities (estimated 
through a pcriodogram) on Figure 2. 

In each case, most of the power lies below /max — 0.005 
which corresponds to a minimum sampling rate of : 

Amin ^ 100 (3) 

This indicates that only 1 parameter on 100 Ci values 
is necessary to reconstruct the whole spectrum. 

Fitting Ct spectra When using the sampling rate A^ = 
Amin, some part of the Fourier spectrum is nonetheless re- 
moved, and the Shannon interpolation given by (2) is no 
longer exact. We therefore modify it to a linear parameteri- 
zation : 



A, 



(4) 



where: 



- the {cij}j=i,...,]Vp,„ represent a set of parameters to be 
estimated on the data through least square minimization. 

- u denotes the "basis" function and will be discussed more 
in details in the next part. 




Figure 2. Power spectral densities (periodogram) of the three Q 
spectra presented on Figure 1. 



- the {Ij} represent a grid oi fixed points: Ij — lo + jAg. 
Note that the first point (lo) can be located anywhere. In 
practice one fixes it to the first measured multipole. 

Eq. (4) denotes a decomposition of the C{1) function 
over the non-orthogonal basis Uj{l) = u[{l — At\. 

Given a set of "data" points, i.e. a set of measure- 
ments and their covariancc matrix, the parameters {sj } and 
their covariance matrix are determined through a simple lin- 
ear least square fit. Then one has an estimate of the C(Z) 
function for any real t value together with its associated 
Gaussian error through standard error propagation (see sec- 
tion 3 for an example). 

Bcises The sampling theorem uses the basis (denoted as 
Shaxmon) : 



u(x) = sine (ttx) 



(5) 



It has the the nice property of being an exact interpolation 
on each point of the grid: 

Ji-li 



u{- 



so that the estimated parameters (oj) satisfy: 

C(lj) = aj 



(6) 



(7) 



and can be estimated "by eye" on a (perfect) C{1) curve. 

However the Shannon basis is not strongly local in real 
space, involving long range correlations in the parameters 
covariance matrix. 

To modify this behavior but still keeping the exact in- 
terpolation property, we multiply it by a Gaussian function 
and call this new basis Gauss-Shannon: 

^2 



u{x) = sine (nx) X cxp ("2^) 



(8) 



With respect to the Shannon basis, the off-diagonal 
terms of the parameters covariance matrix will be reduced. 
Although in principle one has the freedom to adjust the a 



Generic description of CMB power spectra 3 



of the Gaussian, it is more interesting to use "large" values, 
since in that case the bulk part of the Fourier spectrum is 
un-filtered : one just cuts off frequencies near /max and the 
Shannon properties are (mostly) preserved. In the following 
we will therefore investigate the (7 = 3 case in Eq. (8). 

For the sake of comparison, we will also investigate a 
very diflferent basis provided by the simple Gaussian func- 
tion^ {Gauss basis): 



w(a;) =exp(— ) 



(9) 



It does not have the exact interpolation property of Eq. (7) 
since it is always positive. Therefore the estimated param- 
eters ajwill now be very different from the C{lj) function 

values, and the off-diagonal terms of the covariance ma- 
trix will be important. The Gaussian function is however 
the most rapidly decaying in both real and Fourier spaces 
(Hardy 1933) and is therefore of great interest. If we want 
to perform a fit to Eq. (4), we do not have the freedom in 
the choice of a. Indeed, since: 

w(-^) AeU{Aef) oc exp(-27r'(AffT)V') 

by changing a wc just define a new effective sampling rate. 
If wc wish to keep the Shannon critical frequency of Eq. (1) 
we are lead to use cr = 1 to avoid aliasing or over-sampling. 

Edge effects The Shannon theorem, on which is based the 
fit, is defined on an infinite support. By restricting the data 
to some measured £ range, a high frequency ringing appears 
near the boundaries. To reduce this phenomenon one can ex- 
tend the grid beyond the data limits. Given the short range 
of the functions basis used in Eq. (4) , only one or two points 
can be added but this is sufficient in most of the cases. Note 
that we do not add any fake data, but just increase the grid 
position and number of parameters to be determined from 
the fit. 

Results To test the precision of the method, we fit the three 
Ce spectra of Figure 1. We use the Gauss-Shannon basis with 
a grid of 34 points (= 100 and 2 points added on the low-side 
I < 2 and 1 on the high side I > 3000). Since the input data 
has no error, their covariance matrix is defined as identity. 
Figure 3 shows the comparison of the parameterization and 
the genuine spectrum. A close-up of the low £ region is also 
displayed The agreement is very satisfactory up to very low £ 
values. The first few bins arc not perfectly described because 
they introduce a high frequency component that is cutoff by 
the sampling rate of Ai = 100. This can be accounted for by 
a second step decomposition of the residuals as in a wavelet 
analysis. For the sake of simplicity and since experimental 
error bars are high in this region, we choose not to correct for 
this effect in the following. The upper part of the spectrum 
is correctly described. 

To compare the interest of using different bases. Fig- 
ure 4 presents the relative agreement between the fitted 
values and the input WMAP best-fit model for the three 
bases: Shannon, Gauss-Shannon,Gauss. For the modified 

^ Note that this approach is different from the one proposed by 
Odman et al. (2003) where all the parameters of the gaussians 
are fitted. 
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Figure 3. Superposition of the three input Ci speetra (dark thin 
hne) from Figure 1, and the parameterizations described in the 
text (thick grey) . The Gauss-Shannon basis has been used, with a 
sampling rate of = 100 and 2 extra-points on the low side and 
one on the high side (therefore 34 parameters have been fitted). 
The insets show a close-up of the low £ region. For the sake of 
visualization, the thick grey curves have been enlarged. 



bases (Gauss-Shannon, Gauss) , the agreement is better than 
1% in the range £ > 30. For low £ values the agreement de- 
creases up to ~ 15% for £ = 2. 

Both Gauss and Gauss-ShcUiuou give very good results, 
Gauss being slightly more precise for the models we use. 
Recall however that the Gauss basis (unlike Gauss-Shannon) 
provides parameters with values far from the grid points 
and highly correlated. Wc will therefore prefer the Gauss- 
Shannon basis in the following. 

Derived results Once the parameters Ojhave been deter- 
mined, including their covariance matrix Va , one can derive 
easily many characteristics of the spectrum. 

Ce values can be estimated for any (real) I value. Using 
again the notation Uj{l) = u[{l — lj)/A(]: 



(10) 



U being the vector of components Uj = Uj{l). 

The peak (and dip) locations can be determined by set- 
ting the derivatives of Eq. (4) to zero and finding (numeri- 
cally) the roots. 

One can compute also binned values. For a given bin- 
ning (6i, &2, .., feAr), assuming a weighting scheme Wi inside 
the bin b: 



leb 



(11) 
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Figure 4. Relative agreement (in %) between tlie input Ci 
{WMAP best fit model) and the function obtained by fitting 34 
parameters on the Shannon (upper plot), Gauss-Shannon (mid- 
dle) and Gauss (bottom) bases 

with the normalization Wt ~ X^igt™?- 

The covariance matrix of the binned estimate is ob- 
tained again through error propagation B = DVaD^ , D 
being a {N x A'par) matrix defined by: 

Dij =l/WiJ2^i^'jil) (12) 
leb 



3 APPLICATION TO WMAP DATA 

We now turn on to real data, by describing the WMAP TT 
angular power spectrum 

The data consist of a sample of N=899 measurements 
for integer £ values in the [2,900] range, and of the weight 
matrix (inverse covariance) W defined as the Fisher matrix 
for the ML estimate. 

Wc begin with a check of our Ci smoothness hypothesis, 
showing the periodogram of the data on Figure 5. Unlike 
Figure 2, this one has (statistical) noise, described by the 
covariance matrix of the data (C). We estimate the mean 
level of noise power in each Fourier mode k by: 

iV-lJV-l 

< \nk\^ >=J2Y1 expb--^(Z - m)k]Cim (13) 

1=0 m=0 

and show it on the Figure too. 

The (noise subtracted) signal does not show high fre- 
quency components and we keep the cut /max ~ 0.005 deter- 
mined previously, loading to the sampling rate Ai = 100. Wc 
then choose the Gauss-Shannon basis and put 10 Ij knots 

^ first year data vlpl version, as available from 
http://lambda.gsfc.nasa.gov/product/map 
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Figure 5. Power spectrum estimated through the periodogram 
of WMAP data. The dashed line corresponds to the estimate of 
the mean level of noise. 



between 2 and 900. We add two extra^points on the low £ 
side and one on the high side. We are therefore left with 
determining 13 parameters from a set of N=899 input data. 
The (linear) least square estimate reads: 

d = (AWAy^A'^Wc 

(14) 

Va = {AWA)-^ 

with 

Ai,=u{^-^) (15) 

and c is the input data values vector. 

There is a subtlety however. The Fisher matrix (or "cur- 
vature" as defined in Verde et al. (2003)) depends actually 
on the true Ct distribution through the cosmic variance. We 
need a means to incorporate it in our fits, in particular for 
the low t range. This is performed through the following iter- 
ative procedure: wc start from the WMAP best Ci estimates 
(as true values) to obtain the weight matrix, and perform 
the fit. We then recompute the weight from the Fisher ma- 
trix using this time the fitted Ci values, and redo the fit. We 
pursue the iteration until the gets stable. 

For this data set, the fit is stable after 3 iterations. 

Figure 6 shows the result of the last fit. 

The per degree of freedom is excellent: 959/886. Fig- 
ure 7 shows the "pull" distribution that we define as the data 
value minus the estimated value divided by the data error. 
Here again the data errors are obtained from the diagonal 
elements of the inverse of the Fisher matrix. The distribu- 
tion is compatible with a normal one. We checked that this 
feature is valid over the whole £ range. These results indi- 
cate that the data are consistent with the assumption of a 
smooth power spectrum. 

As is clear from Figure 6 our fitted function (which 
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Figure 6. Result of the parameterization performed on WMAP 
data using the Gauss-Shannon basis. Light grey points corre- 
sponds to WMAP data for each integer £. Their error is computed 

from the diagonal elements of the inverse of the Fisher matrix as 
obtained in the last step of the iterative procedure described in 
the text. The dark line corresponds to our best estimate for each 
integer £ and the errors (dark grey area) are computed from the 
parameters covariance matrix through standard error propaga- 
tion. A close-up of the low £ region is also shown. 
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Table 1. Value of the parameters (and grid spacing) obtained 
from the fit to WMAP data on the Gauss-Shannon basis. 



Figure 7. "Pull" distribution of WMAP data w.r.tto our best fit 
estimate. Also superimposed is the fit to a gaussian function. The 
parameters of the histogram and the result of the fit are given in 
the box. 
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Figure 8. Correlation matrix of the parameters determined from 
the fit of WMAP data on the Gauss-Shannon basis. 



eters matrix have also been reduced w.r.t a simple Shannon 
based fit. 

The three bases give however similar results for the final 
fitted function. 

We emphasize once more that this description is inde- 
pendent from any cosmological model. 
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includes the effect of the cosmic variance) has a positive 
slope in the low £ region. 

Finally we provide the parameterization obtained in Ta- 
ble 1 (parameters) and Figure 8 (correlation matrix) . It rep- 
resents the most probable estimate of the band limited Ce 
spectrum and one can derive several secondary results as 
explained in section 2. Note how the parameter values can 
be determined directly from Figure 6, as expected for the 
Gauss-Shannon basis. The off-diagonal terms of the param- 



4 OTHER POWER SPECTRA 

4.1 Fourier transform of rings TmiO) 

For a circular scanning strategy, the Fourier decomposition 
of circles is an interesting analysis method , since it does not 
use map projection and allows to follow in time the detector 
response, avoiding hypotheses such as noise stationarity. 
The power spectrum is simply related to the underlying 
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Ci through (Dclabrouillc et al. 1998) 



(16) 



where 



(i) 9 is the opening angle on the sky (colatitude) 

(ii) Vf^ denote the (normalized) associated Legendre 
polynomials 

(iii) Bi is the beam transfer function, depending just on 
t for a Gaussian symmetric beam. 

The inverse process of retrieving Ce values from a set of 
measured Vm(0) is delicate due to the fact that P™(cos^) 
when t > m/sinO. Ansari et al. (2003) have shown for- 
mally how a Fourier scaling in the flat sky limit allows 
to invert the problem up to very low £ values > 3 for 
{9 = 40°)). The key point in their approach is to interpolate 
the Tm{9) spectrum to non- integer rn values {rn/ sin 9) which 
can be done with the present method, since the mr^(0) 
spectrum is smooth. It is simpler however to stay in the £ 
space and use the previous results. 

We describe again the C{1) function in terms of eq. (4). 
By combining with Eq. (16) : 



E 



2ttB^ , 2 ,£~U 



pr(cos0) w(- 



(17) 



which is still linear in the aj . 

The aj parameters can therefore be fitted directly on 
the Vm (9) data using the above new basis, and the Ce spec- 
trum is still obtained from Eq.(4). 

Note that this approach allows to combine any num- 
ber of detectors with different opening angles and transfer 
functions. 



4.2 Angular power spectrum C{9) 

Using again the decomposition (Eq. (4)) , the two-point an- 
gular correlation function can be expressed linearly in terms 
of the a,-: 



C{9) = Y,aj 



(18) 



The same approach than in the previous part is therefore 
possible. Although equivalent to the usual methods, it allows 
to drop the aim computations and adjust directly the Ce 
distribution, even on partial or masked sky data. 



5 CONCLUSION 

In this article, we take advantage of the smoothness of the 
Ce power spectra to decompose the signal in Fourier space 
and apply an improved version of the sampling theorem. We 
obtain an accurate parameterization of the spectra, that is 
independent from cosmological models, as a function C(Z) 
for real I, with few parameters: for "reasonable" CMB mod- 
els we find that sampling the Ce spectrum with one point 
on 100 is sufficient to retrieve the signal on any £ range at 
the percent level, down to very low £ values. We also show 
how this kind of description can be applied to other CMB 



power spectra as the Fourier spectrum of rings (rm(^)) or 
the angular correlation function {C(9)). 

We apply the method to WMAP first year data and pro- 
vide the complete parameterization of the Ce spectrum. We 
show that the data is consistent with the expected smooth- 
ness of the spectrum. 

Such a parameterization is richer than peaks determinar- 
tion (which can obviously be derived from it) . It also allows 
the combination of various experiments and a test of the 
compatibility between them based on the assumed smooth- 
ness of the spectrum. 

Finally it can allow to compress the data (by a factor 
~ 100) for the storage of cosmological models used in large 
databases. 

The method can be applied to any band-limited func- 
tion to be adjusted on data. 
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